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Classes of spontaneous symmetry breaking at zero and 
low magnetic fields in single quantum dots (QD's) and quan- 
tum dot molecules (QDM's) are discussed in relation to the 
ratio Rw between the interelectron Coulomb repulsion and 
the harmonic confinement, using spin-and-Space unrestricted 
Hartree-Fock calculations. These include: Wigner crystalliza- 
tion for Rw > 1, and formation of non-crystallized electron 
puddles localized on the individual dots in QDM's, as well as 
spin-density waves in single QD's, for Rw < 1. 
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Two-dimensional (2D) electron gases have provided 
(e.g., the fractional quantum Hall effect ||l|,^), and con- 
tinue to provide (e.g., a charge-density wave at higher 
Landau levels [^) a source of discovery of remarkable 
many-body phenomena. Recently, 2D artificial quantum 
dots (QD's) and quantum dot molecules (QDM's) have 
become available, with the capability of controlling the 
dots' size, shape, and number N of electrons 

Single QD's are commonly referred to as "artificial 
atoms", since interpretations of transport and capaci- 
tance experiments draw often on analogies between such 
artificial structures and natural atoms Underlying 
these analogies is an effective (circular) central mean field 
(CMF) picture, with the electronic spectra exhibiting (at 
zero magnetic field) shell closures and following Hund's 
rules for open shells. Indeed, in experiments on single 
QD's, the addition energy (AE) spectra [Q exhibit max- 
ima at the expected closed shells {N ~ 2, 6, 12), and 
at the mid-shells {N ~ 4, 9, and 16) in accordance with 
Hund's rule. 

Here, using the self-consistent spin-and-Space unre- 
stricted Hartree-Fock (sS-UHF) method, we dis- 
cuss, for zero and low magnetic fields (B), three types 
of spontaneous symmetry breakings (SB) in circular sin- 
gle QD's and in lateral QDM's (i.e., formation of ground 
states of lower symmetry than that of the confining po- 
tentials [p2[ ). These include: (I) Wigner crystallization 
(WC) flj in both QD's and QDM's, i.e., (spatial) lo- 
calization of individual electrons, (II) formation of elec- 
tron puddles (EP's) in QDM's, that is localization of the 
electrons on each of the individual dots comprising the 
QDM, but without crystallization within each dot, and 
(III) pure spin-density waves (SDW's) which are not ac- 
companied by spatial localization of the electrons [^). 
Furthermore, we show that CMF descriptions at zero 
and low magnetic fields may apply only for low values 



of the parameter Rw = Q/Tt-loq, where Q is the Coulomb 
interaction strength and Tiloq is the parabolic confine- 
ment; Q = e^/nlo, with k being the dielectric constant, 
^0 = {ti/ 171*100)^^"^ the spatial extension of the lowest 
state's wave function in the parabolic confinement, and 
m* the effective electron mass. With the sS-UHF, we find 
that WC occurs (SB of type I) in both QD's and QDM's 
for Rw > 1. For QDM's with Rw < 1, WC does not de- 
velop and instead EP's may form (SB of type II). We note 
here that while certain quantum-mechanical studies of 
electron localization (WC at high B) in single QD's have 
been discussed previously |p|,|l4|, this is the first study 
to explore, using a self-consistent quantum-mechanical 
treatment, broken symmetry WC and EP states in the 
B = and small B regimes for both circular single QD's 
and lateral QDM's, thus providing new insights into the 
nature of these systems. Additionally, for single QD's 
with Rw < 1 and N < 20, we find that in the majority 
of cases the ground states exhibit CMF behavior without 
symmetry breaking; however, in several instances (e.g., 
iV = 14), a pure SDW (SB of type HI) develops. 

The many-body hamiltonian for our problem is Tt = 
J2iLi -^(*) + Si<j e^/^lri — I , where the second term 
corresponds to the interelectron Coulomb repulsion. The 
single-particle hamiltonian H{i) — -ffo(*) + Kiock(*) + 
Hsii) contains a term describing the motion of an elec- 
tron in a 2D two-center-oscillator (TCO) [|l5| confine- 
ment, i.e., Hoii) = pf/2m* + m*u!Q{xf + where 
y't.i = Vi -yu t = I for t/i < (left) and t = r for yi > 
(right), with yi < and yr > being the centers of the 
left and right oscillator. In the QDM, the dots are joint 
smoothly via a neck described by the term Kicck(*) (4th- 
order polynomial in ^ [[liyi5| ]) which allows variations 
of the bare interdot barrier height (Vf,) for any distance 
d — yr — yi] for d = 0, the system reduces to a single 
harmonically confined QD. Magnetic-field effects are in- 
cluded in i/B(i) = [(p,-eA/c)2-p2]/2m*-h5>BB-S,/;i, 
where A.^ = B{—yi/2,Xi/2,0) and the last term is the 
Zeeman interaction with an effective factor g*, Si is the 
electron spin, and hb the Bohr magneton. To solve the 
sS-UHF equations [|l6|, we use a (variable with d) basis 
consisting of the eigenstates of Ho{i) (due to it's sep- 
arability, the cigcnfunctions of i?o(*) can be expressed 
as products of ID harmonic-oscillator wave functions 
in Xi and parabolic cylinder functions in y'^ ^ [p^). In 
all calculations, we used m* — O.OGTmg (GaAs) and 
hujQ — 5 meV. For k, we used the GaAs value of 12.9 
(i.e., Rw = 1-48), as well as k = 20 (i.e., Rw = 0.95, cor- 
responding to a weakened interelectronic repulsion, due 
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e.g. to the effect of the finite thickness of the dots). 

The spatial distributions of the electronic densities for 
QD's and QDM's reveal in almost all cases that for 
Rw = 1-48 (i.e., k = 12.9) the ground-state solutions 
are Wigner crystallized. An example of such a (finite) 
Wigner crystal with i? = is shown in Fig. 1(a) for a 
d — 70 nm {—yi = ~ 35 nm), Vf, = 10 meV closed- 
sheU QDM with TV = 12. The WC is portrayed by 6 
well-resolved humps (3 in each well) for both the up (t) 
and down (j) spins, and by 6 humps and 6 troughs for 
the spin density (| — J,); note that the density peaks 
for the two spin directions do not overlap. Formation 
of such "Wigner supermolecules" (WSM's) in QDM's is 
analogous to that of "Wigner molecules" (WM's) [Q in 
single QD's. The appearance of such a WC is a con- 
sequence of the large value of Rw (1.48 in this case). 
The mean distance between neighboring density max- 
ima inside each of the coupled dots equals f w 20 nm, 
i.e., roughly twice larger than the effective Bohr radius 
agiK = 12.9) = h'^K/m*e'^ = 10.188 nm. Inspection 
of the wave functions shows that this case corresponds 
to an intermediate electron-density regime, where spa- 
tial localization of individual electrons emerges, but with 
finite-amplitude contributions of each of the wave func- 
tions to several of the density peaks (i.e., "weak" WM, 
see below); full localization into a "classical" WC requires 
even lower densities 

A magnetic field compresses the electronic orbitals in 
the QDM and the consequent increase in Coulomb repul- 
sion promotes electrons to higher orbitals of larger spatial 
extension, with an increase in the spin polarization (spin 
flip) resulting in optimization of exchange-energy gain 
(for a description of such a process in single QD's, see 
Ref. 1^). An example is shown for the QDM in Fig. 1(b) 
for B — 3 T {g* = —0.44), where two of the down-spin 
electrons flipped, resulting in 8 up-spin and 4 down-spin 
electrons, accompanied by a reduced Wigner crystallinity 
(partial "melting") of the WSM, portrayed by the less 
pronounced density peaks [compare Fig. 1(b) with Fig. 
1(a)], and increased density in the interdot region. 

Having discussed formation of a WSM (in a 12e closed- 
shell QDM) made of Wigner molecules in each of the 
coupled dots, we display in Fig. 2 results, with B = 0, for 
the ground state (singlet) of a closed-shell {N = 6) single 
QD for two values of k. For k — 12.9 (i.e., Rw = 1.48), 
we observe again the emergence of a WM; note in Fig. 
2(a, bottom) the six charge-density maxima arranged on 
a ring, with f 20 nm. On the other hand, for the 
same single QD, but with a reduced Coulomb repulsion 
(k = 20, Rw — 0.95), no WC occurs; compare the charge 
densities in Fig. 2(b) and 2(a). Thus by varying Rw, one 
may cross the "phase boundary" separating the localized 
Wigner-crystallization and delocalized (CMF) regimes. 
Furthermore, in the WC regime the electron (charge)- 
localization is accompanied here by a SDW [see top panel 
in Fig. 2(a) and also in Fig. 1(a)]. 



The WM emerging for Rw = 1-48 [Fig. 2(a)] is a 
weak one. The transition into the "strong" -WM regime, 
caused by an increase in the strength of the Coulomb re- 
pulsion, is illustrated in Fig. 2(c) for the same QD but 
with K ~ 6, i.e., Rw = 3.18. The two WM isomers shown 
exhibit sharper electron density peaks reflecting stronger 
localization (seen also from much reduced wave-function 
amplitudes at neighboring sites). The geometries of the 
lower-energy fully spin-polarized, P = N —N j= 6, 
isomer [Fig. 2(c), right], and that of the higher-energy 
isomer with P = [Fig. 2(c), left] agree well with those 
determined classically (i.e., for Rw oo) 

Consider next an open-shell (TV = 6) QDM with d — 70 
nm, Vb = 10 meV, and Rw = 1-48, for which CMF-type 
treatments (as well as local-spin-density functional, LSD 
[prf), predict a total net spin polarization P = 2 in ac- 
cordance with Hund's rule, while we find here that the 
ground-state of the QDM is a Wigner-crystallized singlet, 
i.e. P = [see Fig. 3(a)], consisting of two spin-polarized 
(triplet) WM's [formed inside the left and right dots; see 
Fig. 3(a), top]. An excited state of the molecule (with a 
0.09 meV higher energy) shown in Fig. 3(b) is also crys- 
tallized but with a net spin polarization P = 2; note 
the different spatial configurations of the ground and ex- 
cited states. Reducing the Rw value for the QDM to 
0.95 (i.e., K = 20) transforms the ground-state of the 
6e QDM from the crystallized state [Fig. 3(a)] into one 
consisting of electron puddles [SB of type II, Fig. 3(c, 
left)]; here each of the EP's (on the left and right dots) is 
spin-polarized with P/ = 1, P^ = —1, and the singlet and 
triplet states of the whole QDM are essentially degen- 
erate. Note that the orbitals on the left and right dots 
[see, e.g., those on the left dot in Fig. 3(c, right)] are those 
expected from a CMF treatment, but with slight (ellipti- 
cal) distortions due to the interdot interaction. Only for 
much lower values of Rw (^ 0.20, i.e., k > 90.0) Fermi 
liquid (delocalized) bahavior is restored. 

Results of sS-UHF calculations at P = for the AE's, 
Ae = E{N + 1) - E{N) - [E{N) - E{N - 1)], where 
E{N) is the A^-electron total energy, are shown for single 
QD's and QDM's in Fig. 4. For the single QD's and 
for a rather wide range of Rw > 1 (e.g., see the curve 
marked k = 12.9), the AE's (corresponding mostly to 
WC states) exhibit maxima at the same values of iV (see 
introductory paragraphs) as in the "normal" (i.e., Rw < 
1 non-crystallized, CMF) regime. However, while in the 
latter the spin-polarizations follow Hund's rule (except 
for N = 14, corresponding to a pure SDW state), those 
in the WC states (k = 12.9) for iV > 8 do not. For the 
QDM, a CMF treatment predicts for sufficiently large 
interdot barriers (e.g., Vf, = 10 mcV) shell closures at 
= 4 and 12 (i.e., twice the single QD values) and 
(Hund's) half-shell maximum spin polarizations (e.g., at 
A^ = 8). While shell-closure features are observed in the 
AE's of the QDM's with k = 12.9 (WC) and k = 20 
[where EP's may form, see, e.g.. Fig. 3(c)] shown in Fig. 
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4 (more pronounced for k = 20), in both cases the spin 
polarizations do not in general follow Hund's rule (e.g., 
iV = 6 and 8). 

In summary, using sS-UHF calculations, we discussed 
three classes of spontaneous SB in QD's and QDM's at 
zero and low magnetic fields, i.e., formation of Wigner 
crystallized molecules and supermolecules for Rw > 1, 
and non-crystallized electron puddles localized on the in- 
dividual dots in QDM's, as well as pure SDW's in sin- 
gle QD's (for Rw < !)■ Further studies of such broken 
symmetries may include: mapping of "phase-boundaries" 
through variations of materials dependent (e.g., dielec- 
tric constant) and externally controlled (e.g., gate volt- 
ages, interdot distances and barrier heights, and mag- 
netic fields) parameters, and probing of excitations and 
spin polarizations p9| . 
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FIG. 1. sS-UHF results for a 12e QDM (d = 70 nm, 
Vb = 10 meV, hujo = 5 meV), without [in (a)] and with a 
magnetic field B = 3 T [in (b)]. For both cases, the bottom 
and middle panels correspond to the up-spin and down-spin 
electron distributions, respectively, and the top ones corre- 
spond to the difference between them (spin density). Lengths 
(x and y axes) in nm, density distributions (vertical axes) in 
10~^ nm^^. x-axes, y-axes, and vertical-axes scales in (b) are 
the same as in (a). 

FIG. 2. Total electronic density distributions [bottom of 

(a) and (b), and panel (c)], and spin density distributions 
[top of (a) and (b)] for a 6e single QD with hujo = 5 meV and 
B = 0; K = 12.9 {Rw = 1.48) in (a), k = 20 {Rw = 0.95) in 

(b) , and k — 6.0 (Rw ~ 3.18) in (c). The spin polarizations 
are P = in (a), (b), and (c, left), and P = 6 for the isomer 
in (c, right) whose energy is lower by 1.72 meV than that of 
the one shown in (c, left). Units as in Fig. 1. 

FIG. 3. (a) Ground-state and (b) an excited-state 
Wigner-crystallized total electronic (chaxge) density distribu- 
tions (bottom panels) and spin-density distributions (top) for 
a 6e QDM with Tiloo = 5 meV, d = 70 nm, Vb = 10 meV, 
K = 12.9, and B = 0. (c, left): total electronic density for 
the same QDM, but with k = 20, illustrating formation of 
(non-crystallized within the dots) EP's. (c, right): contour 
plots of the densities of the three individual orbitals localized 
on the left dot {Pi = 1, with spin polarization of the orbitals 
as indicated). Units as in Fig. 1. 



FIG. 4. sS-UHF results for the addition energies {Ae vs. 
N) of a single QD {Tiujo = 5 meV, two upper curves) and for 
a QDM {hcjo = 5 meV, d = 70 nm, Vb = 10 meV) calculated 
for K = 12.9 and 20, and for B = 0. Energies in meV. The 
spin polarizations, P = N —N J,, are marked on the curves. 



4 



Erratum: Spontaneous symmetry breaking in single and molecular quantum dots 

[Phys. Rev. Lett. 82, 5325 (1999)] 

Constantine Yannouleas and Uzi Landman 
(Phys. Rev. Lett. 85, 2220 [4 September 2000]) 



In our recent Letter on spontaneous symmetry breaking in quantum dots (QD's), we displayed in Fig. 4 addition 
energies Ae{N) — E{N + 1) — 2E{N) + E{N — 1), where E{N) is the iV-lectron ground-state total energy, calculated 
with the spin-and-space unrestricted Hartree-Fock (sS-UHF) method. Subsequent to the the publication of our Letter, 
we have performed further systematic sS-UHF calculations, with the use of larger harmonic-oscillator bases and an 
implementation of an extensive search for energy minima. While the behavior and magnitudes of Ae shown in Fig. 
4 of our Letter are maintained, as well as our finding pertaining to the prevalent violation of Hund's first rule, our 
improved calculations yield in certain instances different spin polarizations P = iV f —N J., where iV | and N J, 
{N t +N 1= N) arc the number of electrons with up and down spins, respectively. 

In Fig. 1 below, we display our new results for Ae in a a single QD with GaAs parameters, i.e., k = 12.9, fiujo = 5 
meV, and m* = 0.067me (corresponding to the top curve in Fig. 4 of our Letter). The results of our calculations 
(solid dots) for the spin polarization given in the inset to Fig. 1 exhibit violation of Hund's rule (open squares) for 
N = 4,8, 9, 14, 15, 16, 18, and 22 (note that this violation appears already for = 4); however, as noted in our Letter, 
the addition energies display maxima at closed shells (i.e., at = 6, 12, and 20), as well as at the mid-shell closures 
(i.e., at iV = 4, 9, and 16). In general, the non-Hund ground-state minima are accompanied by energetically close spin 
isomers obeying Hund's rule, and vice versa (e.g., for = 17 the energy difference between the Hund, P = 3, and 
the non-Hund, P = 1, isomers is 0.05 meV). 
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FIG. 5. sS-UHF results for the addition energies (Ae vs TV) of a single QD {hcoo = 5 meV, k, — 12.9, m* = 0.067me). The 
inset displays the calculated (solid dots) spin polarizations, P = N ^ —N I, as well as those (open squares) expected from 
Hund's first rule. 

For the corresponding case of a quantum dot molecule (QDM) with k — 12.9, d — 70 nm, and 14 = 10 meV (third 
curve from the top in Fig. 4 of our Letter), our improved calculation yields P = for A^ ~ 14 (imlike the earlier value 
of P = 2), with all the other spin polarizations remaining unchanged. 
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